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Q\ • Abstract 



Gluon and quark contributions to the thermodynamic potential (free energy) of a 
(2+l)-dimensional QCD model at finite temperature in the background of a constant 
homogeneous chromomagnetic field H combined with Aq condensate is calculated. 
0\ . The role of the tachyonic mode in the gluon energy spectrum is discussed. A possi- 

. bility of the free energy global minimum generation at nonzero values of H and A 

', condensates is investigated. 

00 

g : Introduction 

Quantum field theories in low dimensional space-time have recently excited a con- 
siderable interest [jl], @, |3j due to a close relation with their (3+l)-dimensional 
Ph! analogues [f|, [5], |6j , as well as to a possibility of explaining certain natural phenom- 
ena with their help. For example, the method of dimensional reduction |7j allows 
one to study the QCD3 +i quark-gluon plasma in the high temperature limit above 
the critical temperature of the confinement-deconfinement phase transition. As 
another example, we mention the two-dimensional model of electrons interacting 
with external electromagnetic fields that was employed in the explanation of the 
quantum Hall effect [§, |9| . 

In the last years, much attention has been paid to investigation of the effective 
potential in the background field both at zero and at high temperatures in the 
one loop approximation [|TD|, |TT], [T2], [T3|, |JJ, |T5| , and also with account of higher 



loop contributions []TJ, |T7|, [T8], [2Df . The assumed ^o-condensate formation in 
QCD (see |nj) leads to important physical consequences, such as spontaneous 



breaking of global gauge symmetry, stabilization of the effective potential, elim- 
ination of the infrared divergence etc. New results for the two-loop gauge field 
effective potential in the presence of an external chromomagnetic field have been 
obtained recently in |T9j. At the same time, a problem of gauge dependence of 
the results arises JTS] in the case of nonzero Aq potential, when higher orders of 
the perturbation expansion are taken into account. 
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In the present paper, we study the one-loop gauge field effective potential 
(thermodynamic potential) in the (2+l)-dimensional QCD model in the pres- 
ence of a Aq potential and a chromomagnetic field at finite temperatures. In the 
first section, some basic formulas that describe quantum field systems at finite 
temperatures are reminded, and the relation between the effective potential and 
the energy spectrum of one-particle excitations is demonstrated. In the second 
section, the gluon field contribution to the free energy density is calculated in the 
(2+l)-dimensional model of QCD. The tachyonic mode in the energy spectrum 
of gluons is demonstrated to lead to a nonanalytic dependence of the thermo- 
dynamic potential on the condensate field strength. A possible way of removing 
this nonanalytic behavior is presented. A series of nontrivial minima of the free 
energy, as well as phases of confinement and deconfinement are shown to exist 
when temperatures are below the critical value. In the third section, the quark 
contribution to free energy density is obtained. Calculations show that, in con- 
trast to the gluon sector, its minimum exists in the absence of the condensate 
fields. 



1 Energy spectrum of excitations and free energy 

Consider a QCD2+1 model with a gluon gauge field A a in the adjoint representa- 
tion of the SU (2) color group. It can be written as the sum A® = A^ + a a ^ where 
A a ^ is the background potential and are the gluon quantum fluctuations. The 
gauge field Lagrangian in the Lorentz gauge in the Euclidean space-time has the 
form 

C g = \(FP 2 + \(Dt a V + TiDYx" , (1) 
where = 5 ab d^ — gf abc A^ is the covariant derivative, x> X are ghost fields, 

— 2 t, — etc — cb 

and (D ) = D ^ . In passing to the Euclidean space-time, we retained the 
subscript to denote those components of vectors that are time-like, so that 
xq = it and v = 0, 1, 2. The following representation of the Clifford algebra 
can be chosen in the 3-dimensional space-time: 

7 ° = a 3 ^ 1 = ia\-f 2 = ia 2 , (2) 

where a 1 , a 2 , a 3 are Pauli matrices. Then the Lagrangian for quarks with Nf 
flavors takes the form 



3=1 



N t r \ 



l^d^-ig—AD+mj 



■ (3) 
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Here A a matrices are the color group generators. The generating functional 



Z[A, j, r},rj\= / da a dxdxd^dil) exp 



(4) 



where C = £ g + £ q is the QCD Lagrangian, can be calculated in the one loop ap- 
proximation. To this end, it is sufficient to retain in C g only the terms quadratic 
in gluon fluctuations a®. Namely, 



r - r(°) + r(°) - -CP 1 \ 2 

^g ~~ ' ^g ' ^g — 4^ ' 



£(2) = _I a « 

^g 2 M 



(£ 2 )%, + 2gF c f abc 



(5) 



a 



where = d^A^ — d^A^-^ gf^A^A^. In this case the path integrals in (H) take 
the Gaussian form and can be easily calculated: 



Z\A\ = exp<^ — d z x(F 



[IV J 



det(-L> V - 2gF fiV f) 



x 



X 



det(-L» )n det 

3=1 



A 



(6) 



Using Z = exp (We) 5 one obtains the effective Euclidean action 



W 



E 





2 + (e%) 



(7) 



We have taken into account here that contributions of the gauge field longitudinal 
components and those of the ghost fields cancel each other out. Summation over 
the quantum numbers r and k of the gluon ef and quark e^ k energy spectra 
should include only physical degrees of freedom. Formally, we will consider our 
field system in a three-dimensional cube with volume L 3 and periodic boundary 
conditions. In this case, the effective potential V is determined as follows: 



V = 



W 



E 



L 3 " 



(8) 



The effective potential V in the one-loop approximation can be represented as 
the sum of two terms: V = + v, where = (F^,) 2 /4 is the energy density 
of the background field, v = v g + v q is the sum of the one-loop contributions 
of gluon and quark fluctuations. The limit of an infinite space-time volume is 
obtained, when L — > oo. Since the explicit expression for the effective potential 
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( |T5|) does not depend on L, this limiting procedure will not be mentioned in what 
follows. 

The temperature T = 1/(3 is introduced by imposing periodicity and an- 
tiperiodicity conditions on the boson and the fermion fields respectively in vari- 
able xq with the period equal to (3. According to the well known prescrip- 
tions, this is achieved by introducing the Matsubara frequencies instead of qo, po 
and making the following substitutions: qo — ► 2izl/(3 — i>c\i\ for bosons and 
Po — ► 27r(/ + \)/(3 — %k\i<i for fermions, where are the gluon and quark 

chemical potentials respectively, x = ±1 determines the particles and antiparti- 
cles charge signs, / G Z (Z is the set of integers). The one-loop effective potential, 
determined in (§), is proportional to the free energy density 

Tln(Z) _ n 
~T 2 "Is - 



V = 



(9) 



Upon the above substitution, the quark and gluon contributions to the one loop 
effective potential take the form 

\ 2 



V = V g + V 



1 



q 



1 



-00 



2/3L 2 



E E 



-oo 



Nr 



(3L 2 



EE E 



r, x=±l 



In 




^2 



l>CpL2 



+ (4: 



(10) 



j—1 l=—oo k, K=±l 

We introduce the background field as a superposition of the constant homoge- 
neous chromomagnetic field H, directed along the third coordinate axis, and the 
potential Ao, both pointing in the third direction of the group space. Namely, 



(ii) 



In order to account for the ^o-condensate, it is sufficient to perform the substi- 
tutions — > gAo and i/12 —> gAo/2. 

In the (2+l)-dimensional space-time, fermions are expressed in terms of two- 
component spinors (see (0)), and thus they have no spin degrees of freedom. 
Moreover, the Lorentz gauge fixing conditions a b Q = , D l a b = lead to linear 
interdependence of the positive- and negative-frequency solutions of the Lagrange 
equations 



(D 2 ) ab 5 L 



*abc 



< = 



(12) 



derived from (|5]). In fact, an arbitrary solution of ([T2j) can be represented as an 
eigenvector-series expansion: 



± 



= ai ± ia 2 = ^ N^f n {x), 



(13) 
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where N n are constant coefficients, f n {x) are eigenvectors of equation flT2p . The 
gauge fixing condition D~a + D + a + = can be shown to lead to the following 
restrictions, imposed on the coefficients (see |2T|): N~ = 2{n + 2)N+ +2 , n = 



0,1,2,...; iV^ = 0. No restrictions are imposed on the tachyonic mode coefficient 
Nq. Thus, by setting chemical potentials fii and \i<i equal to zero, we obtain the 
energy spectra of quark and gluon one-particle excitations in the chromomagnetic 
field 

e\ = 2gH(n - i) - ie, n = 0,2,3,4, ... ; (14) 

£q = gHn + m 2 — ze, n = 0, 1, 2, 3, . . . , e > 0. (15) 

Infinitely small negative imaginary term — ie gives a prescription for handling the 
poles, as well as allows to determine a correct limit for the tachyonic mode at 
T -» (see fl22j)). 



2 Gluon contribution to the free energy density 



In order to obtain an explicit expression for the effective potential, one needs to 
substitute ( |T41 ) and ([15]) in f|Kif), with account for the Landau levels degeneracy in 
a homogeneous magnetic field. In the case of the (2+l)-dimensional space-time 
this degeneracy is equal to gHL 2 / (2ir). Then the contribution of the charged 
gluon loop v g to the thermodynamic potential v can be represented in the form 



v b = 



gH_ 

27T/3 



E ln 



l=— oo 



2tx\ 



+ gA ) 2 - gH - ie 



+ 



n=2 



— + gA ) 2 + 2gH(n - -) - ie 



(16) 



where the first term corresponds to the contribution of the tachyonic mode (whose 
energy squared is negative). Limits of summation over / are infinite, and hence, 
expression ( jT^ ) is periodic in gA$ with the period equal to 2n/ (3. This fact is in 
close relation with the gauge invariance property. When T = 0, one can perform 
the following transformation: 

A' = UA Q U + + -Ud Q U + , A\ = UA,U + , where U = exp ^-igx A a ^J . (17) 

Then A' = 0, and so Aq is not a physical parameter. When T^0, the boundary 
conditions A^(xq,x.) = A^Xq + /3,x) lead to subsidiary conditions [t/, A a ] = 
for all A a . By definition, it means that U lies in the center of the gauge group. 
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Therefore, in the case of the SU(2) group, only those gauge transformations are 
allowed that preserve the Z2 symmetry: 

A -> A' = A Q + — — , where n G Z. (18) 



The domain of gauge nonequivalent values of potential Aq is given by gAo G 
[0,2ttT). 

It follows from that v s is real, when the argument of the first logarithm 
in is positive, i.e. under the condition 

< 9^ < 2irT - ^g~H, (19) 

otherwise the vacuum is unstable. As calculations demonstrate, the one-loop 
effective potential in the (3+l)-case has a finite nontrivial minimum, which how- 
ever proves to be unstable [TTj- According to [jTSj, consideration for higher 



loops (ring diagrams) also leads to appearance of an unstable minimum in the 
limit of high temperatures for fields of the order of (gH) 1 ^ 2 ~ g A '^T, which in 
the case of small a s exceeds the one-loop estimates for the field. Returning to 
our general formula (|TB|), we see, that under the condition gA = y/gH, the 



free energy v 8 becomes negative infinite. Such singularity does not arise in the 
(3+l)-dimensional space, because of the smoothing effect of integration over the 
momentum third component, absent in the (2+l)-dimensional case. 

In our case, the divergence can be eliminated by accounting for radiative 
corrections, i.e. due to a nonvanishing imaginary part of the gluon polarization 
operator (PO). We are interested in the energy radiative shift for the tachyonic 
mode only, as it is responsible for the singular behavior of the effective potential. 
In order to make qualitative estimates without any detailed calculations, we write 
down the gluon PO in the form: II(e, T) = a s e(IIi + HI2), where III and II2 are 
some functions of the field and temperature, whose explicit form has no essential 
influence on the further considerations and qualitative results. Note, that II2 
accounts for the gluon decay. Contribution of the PO is of principle significance 
only in the vicinity of the effective potential singularity. Consider the dispersion 
equation s 2 = e 2 + (Il(e,T)}, and make a reasonable assumption, justified by 
further numerical calculations, that the behavior of the free energy near the 
singular point only weakly depends on the values of the functions III ~ II2 ~ 1. 
Then, the energy squared of the tachyonic mode, for a s < (gH) 1 / 2 , with account 
for the radiative corrections and for the equality e 2 = —gH, can be approximately 
found to be: 

<4ch - -9H - iCa s y/g~H, (20) 
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where C is a coefficient of the order unity. As it is seen, the nonvanishing 
imaginary part of £ 2 &ch guarantees the nonzero argument of the first logarithm 
in ( |Tfi| ) at gH ^ 0, though its imaginary part is also nonvanishing. 
Applying the known identity [22] 



n 



1 + 



X 



2irl — a 



ch(x) — cos(a) 
1 — cos(a) 



to (jnj) and discarding an irrelevant additive constant, similarly to [jTT 
obtains 



(21) 



one 



v 



tach 



1 



gH f 
gH_ 

2irf3 



1 



luq + In [1 + e~ 2/3a;o -2 e"^ cos{(3 gA )] j = 



In < 2 cos 



/3(^g~H + -Ca s 



(22) 



2cos(/3gA ( 



V 



tach 



+ i^E{-»V n[1+e 



-2e 



COS 



where 



1 



cj = -iy/g~H + ^Ca s , = 2gH(n - |) 



(23) 



(24) 



An expression for v g can be received with the help of the Fock-Schwinger 
method. It is useful in studying the effective potential in the region of small gH. 
As the first step, recall the standard integral representation 



ln^4 = — f — exp(— sA), 
J s 



(25) 



valid up to an additive infinite constant, and apply it to the second term in QI5Q: 



v 



tach 



gH_ 

2ttP 



OO +O0 p , 

E X — r / (IS 
2_. j y ex P 

n=2 l= — oo n 



— + gA ) 2 + 2gH(n - - 



(26) 



Then, let us separate the function of temperature Vj, in the effective potential 
v g = Vj, = q + Vj, (the zero temperature part Vj, =0 is independent of Aq and will be 
considered below, see Q35f)). To this end, we transform the summation over the 
Matsubara frequencies in (|2TJ) with the help of the following identity (see [[TO 
or 



-co 



^2 ex p 



1= — OO 



,2ttI 
-s[— + gAo)' 



P 



2 x /¥s 



^2 ex p 



2/2 



0*1 



1= — OO 



4s 



■)cos(PgA l). (27) 
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After performing summation over n in fl2§) we obtain the final result: 



.tach 



DO 



ds 



27r 3/2 J s 3/2 




2sh(sgH) 



^exp( 



f3 2 l 



2/2 



1=1 



4s 



cos(j3gA l). (28) 



Numerical estimations confirm that the above expression ( pZ8Q coincides with the 
temperature part of ([23]). The advantage of expression (|28|) is that it has the 
evident limit at gH — * 0: 



v 



T,H=0 



7T(3 



1 X 



1=1 



cos([3gA l) 
/ 3 



(29) 



In the high temperature limit, T ^> gAo, the following approximate expression, 
demonstrating the nonanalytical dependence of the effective potential on the 
background field, can be obtained: 



1 



v. 



T,H=0 



7T/3 £ 



C(3) + 



HPgA Q ) - - 



(30) 



By setting Aq = in (|29|) one obtains twice the effective potential of the 
uncharged gluons (v g contains the contribution of gluons of two opposite color 
charges) 

C(3)T 3 



2v eo = 



TT 



(31) 



As has to be expected, we have obtained an analogue of the Planck law for the 
black body radiation in the (2+l)-dimensional space-time. The total free energy 
density can be obtained by adding v g0 to (|23|) . However, v g0 is the function of 
temperature only, and is independent of the condensate field, and, hence, it is of 
no interest for us. 

In order to obtain numerical estimates of the results, it is convenient to pass 
to dimensionless variables 



x 



f3^/g~H, y = [3gA . 



(32) 



At T = 0, this substitution is evidently incorrect. Nevertheless, in order to 
make our notations universal, we will employ x and y even at zero temperature, 
assuming that in this case (3 takes on a certain finite numerical value in the 
definition (|32|). Let us add a constant to v such that the dimensionless effective 
potential u vanishes when the condensate field goes to zero: 



u(x,y,T) 



v(H,Aq,T) 

J-3 



u(0,0) = 0. 



(33) 



s 



We have to find the minimum of the real part of the function 



x~* T 



U(x, y, T) = U<®(x, T) + u(x, y, T), U®(x, T) = -- 

2 g l 



(34) 



where the quantity T j g 2 provides the scale of the temperature with respect to the 
coupling constant g. Along with the notations u g and u q for the gluon and quark 
contributions to the dimensionless effective potential u, let us also introduce the 
following one: zero-temperature part ut=o of potential u is 

separated, so that u = ut=o + ut and Uj,_ = + Uj, =Q . 

Assuming C = 1, the zero temperature contribution {(3 — > oo) can be obtained 
from 



V T=0 ~ 



gH ^ 



n=0,2 



2tt 



V2 



4tt " C( 2^ 



{gHf 2 , gH 

% V T— Q!s- [dO) 



2tt 



47T 



The real part of this expression coincides with the result obtained in 1)21 



,see 



also The contribution, provided by the gluon tachyonic mode in the ef- 



fective potential, was not obtained in 



Global minimum of U g (Fig. |T]) is 



achieved at y/gH « 0.185^ 2 . There is also a local minimum of Uj< =Q at gH = 
(in contrast to [|2T|1), which is formed by the tachyonic mode contribution pro- 
portional to gH and dominating at gH — ► 0. 

The author of [^T| claims, that the condensate evaporates, when the tem- 
perature is above some critical value T > T CT . The gluon tachyonic mode was 
considered unphysical and omitted in that paper. However, v tach gives nontrivial 
contributions both to the imaginary and real parts of the free energy density at 
nonzero temperature, and hence, it has to be taken into account. Moreover, it is 
the tachyonic mode that can generate a minimum of v g at Aq ^ [pf . 

Examination of the real part of U g (x,y = const) as a function of x reveals 
a nontrivial minimum at that exists at the temperature T lower than 

a certain critical value T cr : T < T CT ~ 0.15g 2 , when the temperature part 
i>f> is small with respect to the zero temperature part wf^o- ^ s an artefact of 
the normalization chosen, the temperature contribution becomes weakly depen- 
dent on temperature, while the zero temperature part becomes a function of 
temperature. In this way, oscillating contribution Uj, (Fig. @) modulates the 
zero-temperature potential U^ =0 , possessing a nontrivial minimum (Fig. [I]). The 
value of x m i n is close to discrete points n7r, n G N and increases with growing tem- 
perature. Moreover, a sequence of second-order phase transitions between phases 
with minima of U g (x m [ n , y) either at y m [ n = 7r or at y m m = takes place. When 
temperature tends to zero, there appears an infinite number of those phases. For 
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example, at T/g 2 = 0.1, a global minimum of the effective potential occurs for 
x = 3.03, y = 7r (Figs. |3|, [|). However, presence of an imaginary part of v g 
means instability of this condensate configuration, at least at the one-loop level 
(cf. When the temperature increases above T cr , the condensate values y m i n 

and x m i n decrease from tt to 0, keeping correlation x m i n ~ y m [ n unchanged. 

As is well known, the system is in the confinement phase, if the trace of the 
Polyakov loop |2j|, [25j] in the fundamental representation vanishes, Trp('P) = 0. 
The Polyakov loop is determined as follows: 



V = T exp 



\a 

dtAl— 



(36) 



In our case ([TT|) potential Aq is directed along the third axis in the color space. 
Hence, Ttf(V) = 2 cos(/5ejo4o/2). It is evident, that the condition Trp('P) = is 
satisfied, when figAo = 7r. Thus, the minima of the effective potential at temper- 
atures below the critical value corresponds to the confinement and deconfinement 
phases. 



3 Quark contribution to the free energy density 



Let us consider the quark contribution to the free energy density in the same 
manner, as it has been done for gluons in section |2j. The quark energy levels 
degeneracy in the external chromomagnetic field is proportional to the quark 
color charge, ±1/2, and is equal to gHL 2 /(47r), which is one half that of the 
gluon case. Upon substitution ([T5|) in ([Itj), the expression for the quark and 
antiquark effective potential is obtained: 



v 1 = 



gH_ 

47T/3 



E E 



i= 



In 

-oo n=0,A=±l 



'7r(2/ + r 



+ gHn + mf 



(37) 



Here A = ±1 corresponds to the quark color charge values. The essential differ- 
ence from the nonabelian gauge field contribution in the chromomagnetic back- 
ground is the absence of the tachyonic mode in the quark energy spectrum. 
Therefore, v q is a well defined function of the condensate field in the whole do- 
main < gH < ex), < gA{) < oo : 



Nf 1 v q = 




n=Q 



where 



1 + e 



u n = gHn 



-2(3u> n 



+2e 



COS 



■PgA 



m 



(38) 



(39) 
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The quark field at finite temperature is subject to the fermionic antiperiodicity 
condition ijj(xo, x) = —i/j(xo + fi, x). Therefore, the period 4ttT of flgg ) in variable 
gAo is twice the period of f[23j) . This implies that the residual gauge 7Li symmetry 



is violated (Z^v in the case of the SU(N) gauge group). In what follows, we 
restrict ourselves to consideration of the quark field in the chiral limit, m = 0. 
In this case the zero temperature quark contribution is equal to 

^o = -g£V^ = ^C(?), (40) 

n=0 

which coincides with the result of [|2T| . In contrast to the case of the gluon 
potential v^^q, the sign of v^ =0 is always positive, and this does not allow a 
nontrivial minimum to be formed. The temperature dependent part of the quark 
effective potential reads 



^4 = -gX> 

M n=0 



x + e -zpu n +2 e -p W „ cog( 



(41) 



By using the proper time method one can obtain an alternative representation, 
which evidently demonstrates continuous behavior of Vj> at gH — > : 

oo 

! q gH f ds ^ (3H 2 cos{(3gA l/2 + irl) 

o 1=1 
In the limit (/iif — > we obtain 

^H^^Ef-!)'^^. (43) 

Z=l 

The dimensionless effective potential u q (x, y), defined in Q33Q , does not depend 
on temperature T. The family of curves for u q with fixed values x = const and 
y = const are depicted in figs. [5|, |6] respectively. It is evident that for any 
x = const the minimum of w q (x, y) is reached at y m [ n = 0. Nevertheless, there is 
a nontrivial minimum at y = const < 3.5 for x m i n ^ 0. Global minimum of the 
gluon effective potential is reached at x = 4.30, y = 0. 

The quark contribution to the free energy is greater than the gluon one at y > 
4, when the quark zero mode plays a dominant role. Therefore, our conclusions 
about the existence of the field condensate and of the confinement-deconfinement 
phase transitions at < y < tt remain unchanged even when the total free energy 
is considered. 
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Conclusions 



Thus, in the present paper, contributions of gluons and quarks to the thermo- 
dynamic potential (free energy) in the (2+l)-dimensional space-time are calcu- 
lated in the one-loop approximation at finite temperature in the background of 
the superposition of the uniform constant chromomagnetic field H and the Aq- 
condensate. Consideration of the role of the tachyonic mode in the gluon energy 
spectrum leads us to a conclusion that it can not be neglected, as opposed to 
what has been done in [jZTj. Moreover, accounting for the one-loop radiative 
correction to the gluon energy spectrum enables a nonanalytic behavior of the 
effective potential, related to the presence of the zero modes in the energy spec- 
trum, to be cured. The free energy minimum is studied and a possibility of its 
formation is demonstrated at certain nonzero values of the chromomagnetic field 
strength H and of the potential Aq. The analysis of the temperature dependence 
of the results demonstrates that below a certain critical value of temperature a 
number of subdomains in the domain of the parameters exist where the system 
is either in the confinement or in the deconfmement phases. This behavior is ex- 
plained by the oscillating contribution of the tachyonic mode to the free energy 
density. Unfortunately, the imaginary part of the effective potential does not 
vanish at points x m [ n ~ 7m, n 6 N, where V(x, y) reaches its minimum. Hence 
the nontrivial minimum of the effective potential generated by the condensate 
fields Aq and H turns out to be unstable. This instability is related to the choice 
of the homogeneous vacuum field, and it is reasonable to assume, that it may 
be removed in the realistic situation of a vacuum field inhomogeneous at large 
distances, where confinement is formed (see also the arguments of [|TTj). With 
regard for what has been said, the states found in the present paper may be 
considered as quasistable. More rigorous justifications for the assumptions made 
above might be obtained in the study of inhomogeneous vacuum field models 
(see, e.g. [pBj), and, following IttSfl, in elaboration of higher loop contributions. 
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Figure 6: Reu q (y). 



